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Abstract. We analyze the direct detection of neutralino dark matter in the 
framework of the Next-to-Minimal Supersymmetric Standard Model. Taking into 
account all the available constraints from LEPII and Tevatron, including bounds from 
the muon anomalous magnetic moment and flavor constraint as b — > sj, we compute 



the ncutralino-nuclcon cross section and compare the results with the sensitivity of 
, dark matter detectors. We also study the relic abundance of neutralinos, comparing 



it with WMAP observations. We show that a light singlet-like Higgs can escape 
accelerator constraints allowing for scenarios with large cross sections, and a very 
interesting phenomenology. 



^ ■ 1. Introduction 



Weakly Interacting Massive Particles (WIMP's) are excellent candidates for dark 
matter, and can in principle be detected via elastic interaction with a nucleus. 
Actually, a large number of experiments are devoted to the direct detection of WIMP's. 
Collaborations like DAMA/LIBRA, EDELWEISS, CDMS, XENON 1 Tonne of Ge/Xe, 
etc. are running and/or in progress pQ. 

Supersymmetric (SUSY) theories with R-parity conservation offer an appealling 
candidate for dark matter: the neutralino as the lightest supersymmetric particle 
(LSP). The theoretical predictions for the direct detection of neutralino dark matter 
in the framework of the Next-to-Minimal Supersymmetric Standard Model (NMSSM) 
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has been analyzed recently [21 [3], and we will use these analyses in what follows. 
Let us remark that the NMSSM provides a solution to the so called ix-problem, 
and makes less severe the " Higgs-little fine tuning" problem in the MSSM, via the 
introduction of a singlet field. Although the symmetries of the NMSSM may give 
rise to the possibility of a cosmological domain wall problem [3J, this can be avoided 
by the introduction of suitable non-renormalisable operators [5]. On the other hand, 
the singlet field mixes with the MSSM Higgses, and the singlino field mixes with the 
MSSM neutralinos, thus offering a richer phenomenology. In particular, it is possible 
to have very light neutralinos. In addition, it is also possible to have a very light Higgs 
which is experimentally viable for a sufficiently high singlet composition. The NMSSM 
superpotential is W = ey (Y u Hi Q l u + Y d H{ Q j d + Y e H\ LP e) - eyA S H{H 3 2 + \nS 3 
and, when the scalar component of S acquires a VEV, an effective interaction \iH\Hi 
is generated, with \i = X(S). We refer to [2] for a detailled analysis of the Higgs scalar 
potential and minimization conditions. 

2. Constraints on the parameter space 

In addition to ensuring the presence of a minimum of the potential, many other 
constraints, both theoretical and experimental, must be imposed on the parameter 
space of the NMSSM: A, k, tan/?, /x, the trilinear soft parameters A x , A K , A UjDjE , the 
gaugino masses M% } 2,3, and the scalar masses mQ^,u,D,E- Let us remark that in our 
computation all these parameters are specified at low scale. A comprehensive analysis 
of the low-energy NMSSM phenomenology can be obtained using the nmhdecay 2.0 
code [B|. We also impose flavor constraints, being the most relevant b —>■ s'j. The recent 
experimental World average for the branching ratio (BR) reported by the Heavy Flavor 
Averaging Group is [7] BR exp (6 — > 57) = (3.55 ±0.27) x 10~ 4 . We include in our code 
the calculation of the b — > 57 branching ratio at the NLO order, following the results 
of [8]. Finally, in our analysis we will also take into account the constraints coming 
from the SUSY contributions to the muon anomalous magnetic moment, = (g^ — 2). 
Note that the current average experimental value is a^ xp = 11659 208(6) x 10~ 10 [9]. 
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3. Results and discussion 

The relevant parameters at low scale are A, k, fi, tan/5, A\, A K the soft gaugino masses 
Mi, M 2 , M 3 . In addition we assume the relation Mi = = ^ that mimics the 
GUT unification. We are interested in the region of the parameter space that enables 
cross section in the sensitivity range of detectors. We further want to be within the 2a 
deviations of the central values of the anomalous magnetic moment of the muon. For 
this reason we choose the relevant slepton parameters Me = 150 GeV, Ml = 150 GeV 
and the trilinear Ae = —2500 GeV. In the quark sector we take Mq d u = Ad = 1 TeV, 
Ajj = 2.4 TeV. We consider that all flavor mixing comes from the CKM matrix. In 
this situation the dominant channel for the b — > 57 is mediated by the charged Higgs. 
The BR(5 — > 37) closely follows the behavior of the charged Higgs mass, which in the 
NMSSM is given at tree level by m 2 H± = f - v 2 A 2 + + . From 

the above, we expect that smaller values of BR(6 — > s 7) should be obtained for large 
k/\ (for positive values of k) or small k/X (for k < 0). In general, smaller values 
of the BR(6 — > s-y) will be also associated to larger values of the product fiA\. On 
the other hand, we have also included in the computation the constraints on the relic 
abundance coming from the combination of the recent WMAP3 observations [10J with 
CMB, large scale structure, and Hubble constant measurements, 0.095 < Qh 2 < 0.112 . 
With respect to the relic abundance, one of the regions where flh 2 is potentially within 
the WMAP range is roughly associated with the band where Mi rs /x [TT]. To fullfil 
the above constraints we also choose Mi = 160 GeV, A\ = 400 GeV, fi = 130 GeV, 
A K = —200 GeV, tan(/3) = 5, as a first example. In Fig{T]we present the BR(6 — > 57) 
(on the left) and the relic abundance (on the right), both in the (A, k) plane. 

We find that for large values of k with respect to A, the LSP is a mixture of 
Bino and Higgsinos. This, together with a comparatively large mass, ensures an 
easy annihilation into gauge boson pairs, so that the relic density is negligible. For 
intermediate values of k, we observe that both the LSP and the lightest Higgs become 
more singlino and singlet respectively. In this example the singlino composition of the 
neutralino is always less than 35 %. As the neutralino mass decreases, some annihilation 
channels become kinematically forbidden, such as annihilation into a pair of Z or W 

t We note that in our computation of the Higgs masses radiative corrections are taken into account [6] . 
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Figure 1. (A, k) parameter space for tan/3 = 5, A\ = 400 GeV, A K = —200 GeV, 
/i = 130 GeV. In both cases the gridded area is excluded by tachyons and the region 
above the thick black line is excluded by the occurrence of a Landau Pole below the 
GUT scale (on the right this region appears in light gray). On the left we plot the 
la and 2cr regions for the BR(& — > sj). On the right we plot the neutralino relic 
density where black points are in agreement with the recents observations mentioned 
in the text, 0.095 < ilh 2 < 0.112 (ID] . For comparison points in agreement with 
old observations, 0.1 < ilh 2 < 0.3 pQ are shown in gray. The vertically ruled area 
corresponds to the occurrence of unphysical minima. 

bosons when m^o < Mz or m^o < My/, respectively. In these cases, flh 2 can be large 
(above the WMAP3 range). Moreover, the spectrum is such that one can encounter 
m h o in 2m^o resonances. Finally, for the regions with very small values of k the Higgs 
becomes lighter than the neutralino and new annihilation channels (the most important 
being XiXi h\h\ and Zh®) are available for the neutralino, thus decreasing its relic 
density. On Fig [2] we plot on the left the neutralino-nucleon cross section as a function 
of the neutralino, and on the right as a function of the lightest Higgs mass. Large cross 
section can be obtained for small Higgs masses m^o w 50 GeV. On Fig. [3] we show 
another case in which it is possible to have a light neutralino with a large component 
of singlino. We display the results for the relic density and the cross section. In this 
case all the (A, k) plane is within the 2a region of accuracy for BR(6 — > 57). On the 
other hand a M is slightly beyond the 2a region. On the left of Fig. [3] we observed a 
region with horizontal lines corresponding to m^o were the neutralino is not the LSP. 
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Figure 2. Scatter plot of the scalar neutralino-nucleon cross section, cXt^o as a 

function of the neutralino mass, m^o (on the left), and as a function of the lightest 
scalar Higgs mass, m h o (on the right). We take A\ — 400 GeV, /x = 130 GeV, 
A K = —200 GeV, and tan/3 = 5. Black (gray) dots correspond to points fulfilling all 
the experimental constraints including the relic density in the range 0.095 < ilh 2 < 
0.112 (0.1 < ilh 2 < 0.3). Light grey dots represent those excluded. On the left the 
sensitivities of present and projected experiments are also depicted with solid and 
dashed lines, respectively. The large (small) area bounded by dotted lines is allowed 
by the DAMA experiment when astrophysical uncertainties are (are not) taken into 
account. 



-4 




Figure 3. For M x = 330 GeV, A\ = 570 GeV, fj, = 160 GeV, A K = -60 GeV, and 
tan/3 = 5. On the left the same conventions as on the right of Fig. 1. In addition 
horizontal lines represent excluded points where the LSP is not the neutralino. On 
the right the same as on the left of Fig. 2. 
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4. Conclusions 

Working in the NMSSM, we have computed the theoretical predictions for the scalar 
neutralino-proton cross section cr^o_ p , and compared it with the sensitivities of present 
and projected dark matter experiments. In the computation we have taken into account 
all available experimental constraint from LEP and Tevatron on the parameter space, 
constraints coming from B and K physics, and the muon anomalous magnetic moment. 
We have also checked the correspondence between the theoretical and observational 
values of the relic density. We have found that large values of cr^o_ p , even within the 
reach of present dark matter detectors, can be obtained in regions of the parameter 
space. This is essentially due to the exchange of very light Higgses. 
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